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^ Abstract 

In previous papers we have presented a calculation describing electrodisintegration of the 
deuteron at GeV energies. The model is fully relativistic and incorporates full spin dependence of 

,__i the final state interactions (FSI), which were obtained from the SAID analysis. It was, however, 

> 

f^ limited kinematically due to lack of availability of the SAID amplitudes. This work rectifies this 

problem by implementing a Regge model to describe the FSI. We present an outline of the model 

o 

t-*^ and show comparisons between the two approaches in a region of overlap. We see good agree- 
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ment between the models, and note observables which can provide additional insight due to model 



sensitivity. 
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I. INTRODUCTION 

In previous papers we introduced a new model of deuteron electrodisintegration for use 
at Q 2 > 1 GeV 2 . The model is based on the Bethe-Salpeter equation and uses the distorted- 
wave impulse approximation (DWIA) with the bound state vertex function used in previous 
calculations of elastic electron scattering from the deuteron and a one-body current of the 
usual Dirac-plus-Pauli form. The inclusion of final state interactions used in this model 
contain all possible spin-dependence of the np scattering matrix unlike previous calculations 
intended for use at large momentum transfers. Calculations were presented for observables 
for unpolarized hadrons PQ, for deuteron polarization [2] and for polarization of the ejected 
proton [3]. 

The final state interactions are included in the model by means of the Fermi invariant 
parameterization 



m = Js(s,t)i (1) i (2) + JHMh M1) 7i 2) + -MMK K1 V 



(2) 
jiv 



-Jp( S ^)(z 75 ) (1) (^5) (2) H-^4( S ,0(757 m ) (1) (757m) (2) (1) 

which includes all possible independent spin structures for onshell NN scattering. In the 
previous work, the invariant functions J-i(s,£) were obtained by calculating the five indepen- 
dent helicity matrix elements of (fTl), inverting these equations to obtain the J-"i(s, t) in terms 
of the helicity matrix elements and using the np helicity matrix elements available from the 
SAID partial wave analysis [U |5]. Since the SAID helicity amplitudes are only reliable up 
to a laboratory kinetic energy of 1.3 GeV, this placed limits on the kinematics that could 
be reached by the calculations. 

In order to eliminate this limitation, we have recently completed a fit to available NN 
cross sections and spin observables from s = 5.4 GeV 2 to s = 4000 GeV 2 using a model 
based on Regge theory j6] . The purpose of this paper is to compare the results for deuteron 
electrodisintegration observables using the final state interactions based on the SAID and 
Regge model versions of the final state interactions in a kinematic region where the two 
methods overlap. 

In the next section we provide an outline of the production of the Fermi invariants 
from the Regge model. The following section contains the results of the calculations of 
the electrodisintegration observables and a discussion of these results for the two final state 
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FIG. 1. Isoscalar (a) and isovector (b) mesons with spin J plotted versus the square of the meson 
masses \x. The various lines correspond to the Regge trajectories used in the fit to NN scattering. 
An additional trajectory, the Pomeron, with ++ is required to fit the large s data. It has an 
intercept of 1.08 and a slope of 0.25 GeV 2 . 

interaction models. The last section contains conclusions based on these calculations. 
II. THE REGGE MODEL 



There is an extensive literature describing the origins and applications of Regge theory. 
We provide a short summary of this approach as used in |6] for those who may not be 
familiar with material. A more extensive introduction can be found in [7HTT]. 

The motivation for Regge theory in the case of NN scattering is the observation that 
mesons with fixed parity, G-parity and isospin when plotted for the spin J versus the square 
of mass appear to lie along smooth curves as shown in Fig. [TJ The curves can be represented 
as 

J = a(/i 2 ) (2) 

where a(/i 2 ) is some function of the square of the meson mass /i 2 . In the case of the well 
established mesons, the function is consistent with a straight line. The interpolating func- 
tions «j describe Regge trajectories. Regge theory describes the NN scattering amplitudes 
in terms of the exchange of Regge trajectories rather than individual mesons. 

For NN scattering the application of Regge theory can be simplified by noting that the 
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(a) 

FIG. 2. (color online) Pictorial representation of the direct and exchange contributions to the NN 
scattering amplitude. 

total scattering amplitude can be described in terms of a direct contribution, Fig. |2](a), and 
an exchange contribution Fig. [21(b), where any set of possible diagrams represented by the 
ellipse is the same in both cases. This is due to the fact that any exchange of nucleons in the 
summation scattering diagrams can always be "unwound" to give contributions involving 
the exchange of nucleons in the final state. For this reason, we can concentrate on Fig. [2Fa) 
since Fig. J2?b) can be obtained by the replacement of Mandelstam t with u and a mixing 
of amplitudes due to the interchange of spins in the final state. The scattering operator 
associated with Fig. [21(a) can be written in Fermi invariant form as 

M(«) = F|(M)1 (1) 1 (2) + i^(M)7 K1) 7£ 2) + Fl(9,t)avW(T$ 

-Fi( S ,t)(i 75 )( 1 )(i 75 )( 2 ) + Fl( S ,t)( 757 ^( 1 )( 757M )( 2 ). (3) 

Now consider a contribution to Fig. [21(a) from the exchange of a single meson. The meson 
has definite values of parity (P), G-parity (G), isospin (I) and spin (J). The propagator for 
the exchanged meson is a function of t, and since in the s-channel center of momentum (cm) 
frame t < 0, the pole in the propagator can not be reached. For the purposes of Reggeization, 
it is useful to consider this contribution in the t-channel cm frame as represented by Fig. [3j 
This represents NN annihilating to produce a meson which then converts back to an NN 



FIG. 3. (color online) Representation of the direct contribution to NN scattering through a single 
meson exchange in ^-channel center of momentum frame. 

state. In this frame t > Am 2 and s < 0. So for mesons of sufficiently large mass, it is possible 
to reach the pole in the propagator. The contribution of this diagram can be expanded in 
the partial wave series in the helicity basis. 



<i™M = £< w + 1) feU*<«> - P(-rr' +x ^'f"iU.^(ft 



j=o 



d\> 2 -\ 2t \> 1 -\ 1 (&t) , 



(4) 

where the Aj are the helicities of the respective nucleons, /y 3 / y Al (s) is the scattering 
amplitude in the helicity basis for a given /, P, G and J. 

We next construct initial and final spinor states of good parity, G-parity and isospin and 
use these to produce matrix elements of the Fermi invariant expression (J3J). Comparison of 
this result shows that only helicity amplitudes with A' 2 = A2 and X[ = Ai are required to 
obtain the Fermi invariants F*(s,t). So 



M 



SLaJM) = $> J + !) [fliMMM^ ~ n-lY f^ X2MM (t)} PM , (5) 
j=o 



where we have used dQ (9 t ) = Pj(z) and 

z = cos9 t 



2s 



1 



(6) 



Am 2 - t 
Furthermore, the Fermi invariants can be written in terms of the particular combinations of 
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FIG. 4. (color online) Integration contours used in Reggeizing the scattering amplitudes. Diagram 
(a) represents the original contour used to replace the original sum over integer J with an integral 
over complex J. The dashed line indicates that the contour is closed at infinity. Diagram (b) shows 
the distorted contour which picks a Regge pole through contour C" . 

helicity amplitudes defined as 



M^ )IPG (s,t) = MJTHM) =F PMTri i 



r (a)IPG 

1111' 
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(a)IPG 

l _i j 
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V(2 J + 1) Ul^i j {t) t PfLFA i i (*)1 [1 ± (-1) J ] Pj{z) 

*— ' L 2'2'2'2 2' 2'2'2 J L J 

J=0 
oo 

^(2J + l)/I^(t)[l±(-l) J ]P^) (7) 



j=o 



Note that all of the dependence on s is contained in the Legendre polynomial through the 
definition of z. 

The scattering matrix can now be Reggeized by replacing the sum in (Is]) by performing 
an integral over the contour C, as shown Fig. |4Fa), to give 



AL 



(a)IPG 



± 



1 f dj (-lY(2J + l) fl PG (J,t) [1 ± (-1)'] Pj(z) 



2/ 



c 



sin7rJ 



I f dJ (2J+l)fi P G (J,t)(P J (-z)±P J (z)) 



1i J c sin7rJ 



(8) 



The contour is next distorted as shown in Fig. W[b). Assuming that there is at least one 
pole in the complex- J plane at J = a(t) introduces the closed contour C" and assuming that 



the semicircular part of C vanishes at infinity gives 



2iJ_i_ ioo simiJ 

| 7r(2a(t) + l)^ PG (t) (P aft) (-z) ± P a(t) (s)) 
sin7ra(t) 

This is the Sommerfeld-Watson transform. The first term is the background integral which 

is unknown. The second term describes the Regge pole contribution where /3± p (t) is the 

residue of f± PG (J,t) and a(t) is the trajectory of the Regge pole. 

This expression can now be analytically continued to the s-channel cm frame which 

describes physical A^A" scattering. In this frame s > 4m 2 , 4m 2 — s < t < and z > 1. For 

z > 0, 

2 sin TTfv 

P a (-z) = Q a (z) + e- l ™P a (z) (10) 

7T 

And for z > > 1 the asymptotic forms of the Legendre functions are given by 

V7rr(a + 1) 

Qa(z) -y AI Y + V (a+1) - (I 2 ) 

For large values of z, the background integral goes as z~z. So this contribution can be 
suppressed by requiring that z be large which implies that s must be large compared to 
4m 2 and \t\ must be small. The Legendre function in the Regge pole contribution can also 
be expanded in z. Since the residue /^ it) is undetermined, we can absorb all of the 
real t dependent factors in the second term of (J9J) and those arising from the expansion 
of the Legendre functions into the residue function /3j_ PG (£). The Regge pole contributions 
necessary to determine the invariants -F 7 (s,t) can be parameterized as 

Rg°{8,t) = C(s,t)J2tk±(t)Pl P k G (t> ah{t \ (13) 

k 

where (3 IPG (t) is the residue and £±(£) is a phase function. In fitting the Regge model to 

A^ A" scattering observables, we use three forms of the residue function given by 

0i(t) = foe** 

^ II (t) = (l-e' t )^ lt (14) 

Pin® = -^M 11 



where /3q, Pi and 7 are fitting parameters. The phase function is 

{-i(na(t)/2+S) , 

(15) 
e -i(7r(a(t)+l)/2+5) _ 

where a(t) = a^ + a±t, and 5 is a phase that is a fitting parameter. The function ((s,t) is 
a cutoff factor defined as 

C(M)= A-e 20 (3^l- 1 )') . (16) 

This was introduced to improve the stability of the fitting procedure near 6 = 0° and 180°. 
It has no effect at large s. 

The five classes of Regge poles needed to produce the invariants F?(s,t) are 

R[(s,t) = R I + + (s,t) 
Ri(s,t) = R I _Y(s,t) 

R I 3 (s,t) = R I+ f(s,t) 
Rl(s,t) = R I + Y(s,t) 
R I 5 (s,t) = R I + - 5 + (s,t) (17) 

and the invariant functions are given by 

F t I (s,t) = ^(M)4(M), (18) 

where i G {S,V,T,P,A} and 

E$i(M) : 
Sy 2 (M) : 



2 
m 



2(4m 2 - 1) 
Am 2 -t 



8(2s + t- Am 2 ) 

St3(s,*) = " 



2 

m 



4(2s + £ - 4m 2 ) 

5 
P4is,£) = -- 



2 



2t 
l A5 (s,t) = l. (19) 



All other components of the matrix are 0. 



We also utilize a factor of -^ for type 4 exchanges, which was necessary to guarantee 
conservation of angular momentum. Type 5 exchanges were multiplied by ^-, which we 
assume we can factor from F^s^t). This is necessary in order to cancel with an additional 
factor of s, which occurs when performing the Dirac algebra to calculate the NN helicity 
amplitudes, and ensures unitarity for large s. 

Parametrizing the Fermi invariants in terms of Regge poles provides a straightforward 
method of ensuring that Regge exchanges of appropriate quantum numbers contributes 
appropriately to the NN scattering. It also greatly simplifies the Regge analysis in that it 
reduces to the spinless case. While the techniques vary, conceptually similar approaches 
have been used in the past [12], [13] . 

The above describes only the contributions of the direct diagram Fig. [21(a). The exchange 
contributions, Fig. [21(b), are obtained by the replacement t — >• u and a mixing of the 
invariants. The complete pn invariants can then be written as 



JT(s,t) = F° k (s,t) - Ffo,t) - 2S u Ffo, 



u 



(20) 



where 
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The complete pp invariants are given by 
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(21) 



JFf{s,t) = F° k (s,t) + Fi(s,t) - S kl (Ff(a,u) + Ffa,u)) 



(22) 



Our object was to construct a fit to A^ A" scattering observables that could be used at higher 
invariant masses than allowed by the SAID partial wave analysis. We were motivated to start 
with Regge theory as a result of its success in describing this process for large values of s and 
then to connect it to the region of 5.4 GeV 2 < s < 12 GeV 2 where a considerable amount of 
data is available for a variety of observables. In doing this, we are clearly venturing into the 
region where many of the assumptions of Regge theory are no longer valid. Our approach 
was to retain the form of the Regge theory to provide fits in this region by introducing a 
set of "effective" trajectories with no relation to the meson trajectories to provide a basis 



TABLE I. The kinematics for the x = 1 and x = 1.3 kinematics used in the electrodisintegration 
calculations presented here. 



kinmatics 


1 2 


X 


1.0 1.3 


Q 2 (GeV) 2 


2.4 4.5 


Ebeam (GeV) 


6.25 8.6 


s (GeV) 2 


5.91 5.93 


e (deg) 


15.97 16.0 


Q 2 


0.60 0.57 



for fitting in this region. The fit which we obtained was respectable given the wide range of 
s and the large number of data sets produced at a variety of institutions over a long time 
period, but may obviously be subject to improvement. We will now compare the final state 
interactions contributions to deuteron electrodisintegration observables for this fit and that 
derived from the SAID helicity amplitudes. 

III. RESULTS 

In this section we present a comparison to various observables for deuteron electrodisin- 
tegration for unpolarized hadrons, polarized deuteron and polarized final state proton. Two 
different kinematics are used, one for Bjorken x of x = 1.0 and the other for x = 1.3. Various 
kinematical variables for the two kinematics are shown in Table [H The kinematics are chosen 
such that s, the electron scattering angle 6 e and the ratio of the square of four-momentum 
transfer to three-momentum transfer Q 2 /q 2 are approximately equal for the two cases. The 
values of s are close to the upper range available from SAID and are at the lower end of 
the fitting range for the Regge parameterization. In all cases the onshell approximation 
for the final state interactions (FSI), as described in [T], is used. An offshell prescription 
for the SAID FSI was proposed in [T], but a more complete approach is possible for the 
Regge parameterization and will be considered in a future paper. Since both the SAID and 
Regge amplitudes are fit to onshell data, using only the onshell approximation is the most 
reasonable way to compare the two methods. 

10 



Figure [5] shows the observables for the case where neither the deuteron target nor the 
ejected proton are polarized, and where the azimuthal angle is chosen to be = 180°. A sum- 
mary of the various calculated quantities is contained in the appendix. Figures |5^a) and (b) 
show the differential cross sections as a function of missing momentum p m for the plane wave 
impulse approximation (PWIA) and for the SAID and Regge FSI for the x — 1 kinematics 
and the x = 1.3 kinematics respectively. The size and shape of the two FSI calculations 
are similar in each case. Since these are semi-log plots, a more accurate evaluation of the 
differences is given by the ratio of distorted wave to PWIA cross section a rat i as is shown 
in Fig. [6] In Fig. [6ta) for the x — 1 kinematics, the SAID and Regge results are very 
similar for p m < 0.3 GeV but differ by up to 50 percent from the PWIA result. At higher 
missing momenta both the SAID and Regge results become increasingly large compared to 
the PWIA, and reach a value of approximately 8 times the PWIA at p m = 1.0 GeV for the 
SAID FSI and approximately 5 times for the Regge FSI. For the x = 1.3 kinematics, shown 
in Fig. [6tb), the difference between the SAID and Regge FSI are much smaller and they are 
both much closer to the PWIA. Note that for x = 1 both of the FSI lie above the PWIA 
but for x = 1.3 they are below. This suggests that it may be possible to find kinematics 
at which the FSI effects are minimal and may allow for an approximate extraction of the 
deuteron ground-state momentum distribution, as has been suggested previously [Hj. 

Figures [51(c) and (d) show the transverse-transverse asymmetry Att for the x = 1 and 
x = 1.3 kinematics, respectively. This asymmetry, which is proportional to the Rtt response 
function, is generally assumed to be small since Rtt has generally been shown to be small. 
This is the case in [5Fc) for the PWIA and SAID results, but the asymmetry for the Regge 
FSI is large for intermediate values of p m . The reason for this can be seen from Fig. [7] which 
shows Rtt for the x = 1 kinematics. Note that all three calculations minimum at around 
p m = 0.05 GeV where the cross section is large. However, while the PWIA and Regge 
results fall smoothly to with increasing p m , the Regge results show a second minimum in a 
region where it is comparable in magnitude to the rapidly falling cross section. This results 
in the large values for Att which involves a ratio of the transverse-transverse contribution 
to the cross section to the sum of the longitudinal and transverse contributions. It should 
be noted that the relationship between the Fermi invariants and the response functions is 
very complicated and can involve interferences between the various contributions. As a 
result we have not been able to isolate a single source for the second peak in the Regge Rtt 

11 
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FIG. 5. (color online) Spin observables for unpolarized hadrons. Short-dashed lines represent the 
PWIA contribution. Long-dashed lines include the SAID FSI and solid lines include the Regge 
FSI. Plots in the left-hand column are for the x = 1 kinematics and plots in the right-hand column 
are for the x = 1.3 kinematics. 
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FIG. 6. (color online) Ratios of the differential cross sections to the PWIA approximation. Lines 
are represented as in Fig. [5} 
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FIG. 7. The transverse-transverse response function Rtt for the x = 1 kinematics. Lines are 
represented as in Fig. [5} 

response function. Interference response functions and their associated asymmetries may 
show unpredictable sensitivities to small differences in the Fermi invariants. Ascertaining 
the significance of these differences requires that the errors in fitting parameters for the 
scattering amplitudes be propagated to the electrodisintegration calculations. This can be 
done for the Regge case since we can generate the hessian matrix for the fit. This will be 
done when sufficient resources are available. Unfortunately, we do not have access to similar 
information about the SAID helicity amplitudes. 

Figures |5te) and (f ) show the longitudinal transverse asymmetry Alt for the x = 1 and 
x = 1.3 kinematics respectively. At x = 1, this asymmetry is relatively large and the two 
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FSI models give comparable results and differ substantially form the PWIA. At x = 1.3, the 
two FSI models have similar form but tend to be in less agreement than in the x = 1 case. 
Both, however, are much closer to the PWIA result. 

Figures [5tg) and (h) show the longitudinal transverse asymmetry Alt 1 for the x = 1 
and x — 1.3 kinematics respectively. Measurement of the asymmetry requires a polarized 
electron beam. Since this response is odd under the combination of time reversal and parity, 
its value is in PWIA. For both kinematics it is small for both FSI models. The significance 
of the differences between the SAID and Regge results is unclear. 

Figure [8] shows the single and double spin asymmetries for vector and tensor polarization 
of the target deuteron along the direction of the electron beam at an azimuthal angle of 
<p = 35°. We see good agreement in these observables between the Regge and SAID approach 
as well as strong effects from the FSI. This suggests that target polarization asymmetries 
can provide insight to the effects of FSI while masking the model dependence of how these 
are calculated. 

Figures p[a) and (b) shows the vector polarized target asymmetry A d , for x = 1 and 
x = 1.3 kinematics respectively. Note that this asymmetry is zero in the absence of final 
state interactions. Qualitatively the Regge and SAID approaches are similar. 

Figures p[c) and (d) show tensor polarized target asymmetry Aj, for x = 1 and x = 1.3 
kinematics respectively. Here the two approaches are in excellent agreement and we see a 
dramatic change in behavior for x — 1. The FSI contributions to this observable are minimal 
at the x = 1.3 kinematics. 

Figures M (e) and (f ) show the double spin asymmetry for vector polarized target and 
polarized beam A^ d , for x = 1 and x = 1.3 kinematics respectively. We again note that 
the two approaches are in excellent agreement and observe that the FSI contributions are 
minimal at the x = 1.3. 

Figures [81(g) and (h) show the double spin asymmetry for tensor polarized target and 
polarized beam Aj d , for x = 1 and x = 1.3 kinematics respectively. This asymmetry is zero 
in the PWIA. Qualitatively the approaches yield similar results, and while the FSI do cause 
a non zero contribution the value is relatively small. 

In Figure [9] we present the results for polarized ejected proton at an azimuthal angle of 
= 35°. All asymmetries in Figure [9] are zero for PWIA, thus presenting an ideal set of 
asymmetries for exploring the contribution of FSI. 
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FIG. 8. (color online) Single and double spin asymmetries for vector and tensor polarizations 
along the beam axis. Plots in the left-hand column are for the x = 1 kinematics and plots in the 
right-hand column are for the x = 1.3 kinematics. Lines are represent as in Fig. [5j 
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Figures 9[a) and (b) show the asymmetry A™ , for x = 1 and x = 1.3 kinematics re- 
spectively, and we see good agreement between the models. Figures |9|c) and (d) show the 
asymmetry A 1 . For the x = 1 kinematics we observe that the model approaches are similar 
in magnitude but differ in sign. Because of the strong model dependence evident in this 
observable, and due to the relatively large value, this would provide an interesting mea- 
surement, which could shed light on the role of FSI as well as the various models used to 
calculate them. FSI effects at x = 1.3 are less pronounced. Figures [91(e) and (f ) show the 
asymmetry A s . Here we again see good qualitative agreement between the two models. 



Figure [10] shows the double spin asymmetries for polarized beam and polarized ejected 
proton at an azimuthal angle of = 35°. Figures 10 a) and (b) show the asymmetry A™ , 



for x — 1 and x — 1.3 kinematics respectively. Again we note that this asymmetry is highly 
sensitive to FSI model dependence at x = 1 causing deviation in opposite directions to the 
PWIA, although the magnitude of the deviation is relatively small. The same behavior is 
observed for x — 1.3, however less dramatic. 



Figures 10 c) and (d) show the asymmetry A 1 , for x — 1 and x = 1.3 kinematics 
respectively. In this case we observe qualitatively the same behavior between the SAID and 
Regge approaches although the Regge model is much more drastic at x = 1. At x = 1.3 
we see similar, albeit less pronounced effects. Due to the large differences between the 
approaches we again point out that measurements of this asymmetry would prove insightful. 



Figures 10 [e) and (f) show the asymmetry A s , for x — 1 and x = 1.3 kinematics 
respectively. We note that for both kinematics the two models are qualitatively similar, 
with relatively small deviations from the PWIA and each other. 

IV. CONCLUSIONS 

In this paper we have presented a new method of calculating final state interactions for 
the electrodisintegration of the deuteron. The FSI are calculated using a Regge model, which 
was fit to available NN scattering data. The model is fully relativistic and incorporates full 
spin dependence. With the addition of this new method we have significantly extended 
the kinematic range of our deuteron electrodisintegration calculation which was previously 
limited due to the absence of high energy NN amplitudes from SAID. We have presented 
results for kinematic regions where the SAID and Regge approaches overlap for comparison. 
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FIG. 9. (color online) Single spin asymmetries for ejected protons polarized along the h' , I' and s' 
directions. Plots in the left-hand column are for the x = 1 kinematics and plots in the right-hand 
column are for the x = 1.3 kinematics. Lines are represent as in Fig. [5j 

The comparisons suggest that for most of the observables there is good agreement be- 
tween the two approaches. We expect that most discrepancies between the two models 
would be within error bands were they available. We anticipate being able to propagate the 
error for the Regge model once sufficient resources are available, however, an error analysis 
requires information that is unavailable from SAID. Propagation of the NN fitting error 
to the electrodisintegration observables will require a substantial amount of computational 
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FIG. 10. (color online) Double spin asymmetries for ejected protons polarized along the h' , V and 
s' directions. Plots in the left-hand column are for the x = 1 kinematics and plots in the right-hand 
column are for the x = 1.3 kinematics. Lines are represent as in Fig. [5j 

resources. 

The results are consistent with expectations that FSI play a vital role in understanding 
the reaction mechanism. We have noted that from our results there may be kinematic 
regions where FSI are minimized and PWIA is a valid approximation. This is most evident 
in the cross section and polarized target asymmetries. We have also identified asymmetries 
with large FSI contribution and significant sensitivity to the model dependence of the two 
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FIG. 11. (Color online) Coordinate systems for the D(e,e'p) reaction, k and kl are the initial 
and final electron four-momenta, q is the four-momentum of the virtual photon and p is the four- 
momentum of the final-state proton. 



approaches. In particular Att, A l p and A l ep . While almost all observables are sensitive to 
FSI and measurements would prove useful, these are particularly interesting because of the 
discrepancies between the two models. 



Now that we have tested our Regge model in the region of overlap available with SAID 
we anticipate future work. Offshell effects should be taken into account and the Regge 
model allows for a natural offshell extrapolation. In addition we now have the capabilities 
to explore the kinematic region which will become accessible once the Jefferson Lab 12 GeV 
upgrade is completed. 



Acknowledgments: This work was supported in part by funds provided by the U.S. 
Department of Energy (DOE) under cooperative research agreement under No. DE-AC05- 
84ER40150 and by the National Science Foundation under grant No. PHY-1002478. 

19 



Appendix A: Deuteron Electrodisintegration Observables 

The hadronic response tensor for deuteron electrodisintegration with polarization of the 
target deuteron and the ejected proton can be written as 

w^Xy(D,S) = Y^ (pis[;p2S 2 ; (-)| Ja; \PX' d )* (pis 1 ;p 2 s 2 ; (-)| J\ 1 \PX d ) 

Px d y d (D)P p Sls[ (S) (Al) 



where 



and 



J ±1 = =f-^( J 1 ± J 2 ) (A2) 



J = J°. (A3) 

The deuteron density matrix p® ,, (D) is a function of 

Pe{C/,T 10 ,T 11; T 2 o,T 21 ,T 22 }, (A4) 

where U designates an unpolarized deuteron and the T^ are the components of the polar- 
ization tensor as described in more detail in [2]. The proton density matrix is given by 

P P S1S ,(S) = 1x1(1 + <? ■ S)X S[ ■ (A5) 

where S is a unit vector in the direction of the proton polariation it its rest frame. 

The general form of the D(e, e'p) cross section can be written in the lab frame as [151 US] 



da 5 



dE'dn K dn„ J u n c 16tt 3 M d 



TflpTfl n Pp j 

m q ,, G Mott Jrec 






vlK l + vt7^t + vttTItt + v lt TIlt 

(A6) 



where M^, m p and m n are the masses of the deuteron, proton and neutron, p p = pi and fi p 

are the momentum and solid angle of the ejected proton, E' is the energy of the detected 

electron and Q e is its solid angle, with h = ±1 for positive and negative electron helicity. 

The Mott cross section is 

_ / «cos(fl e /2) \ 2 

\2esm 2 (6 e /2)J 
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and the recoil factor is given by 



fr 



upp — E p q cos 6 P 



M d p p 



(Ah 



The leptonic coefficients vk are 



v L 



Q 2 



r/:= 2? +tan2 y 



vtt = — 



2g 2 



Q 2 Q 2 Mi 2 o e 



VLT' 



Q 2 + Oe 

—pz — tan — 
vV 2 



6 e Q 9 6 e 

v T > = tan — a -— + tan z — 

2 y g 2 



(A9) 
(MO) 
(AH) 

(A12) 
(A13) 

(A14) 



The response functions are related to the hadronic tensor by 



K L (D,S)= Woo (D,S) 

TZ T (D, S) = w n (D, S) + w_i_i(D, S) 
Ktt(D,S)=2$1(w 1 - 1 (D,S)) 
K LT (D,S) = -2&(w 01 (D,S)-wo-i(D,S)) 
n LT ,(D, S) = -2U(w 01 (D, S) + w„-i(A S)) 

n T ,(D,S) = Wll (D,S) - w_i_i(£>,5) . 



(A15) 
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1. Unpolarized 

In the case where the deuteron and proton are not polarized the response functions can 
be written as 

n L (u,6) = R L 

Kt{U, 6) = R T 
n TT (U,0) =R TT cos2(f> 

n LT (u, 6) = r lt cos <\> 

K LT , (U, 6) = R LT > sin 
Kt>(U,0) = (A16) 

where the dependence on the azimuthal angle is written explicitly and the response functions 
Ri are independent of this angle. 

Three asymmetries can be defined for the interference response functions, 

Att= ^tRtt 

A LT = ^^ (A18) 

vlRl + vtRt + vttRtt 

and 

vlRl -\- vj'Rx — VttRtt 
Note that while Alt can be obtained by measuring protons in the electron scattering plane 

symmetrically about the direction of the three-momentum transfer, the asymmetries Att 

and Alt 1 require measurements to be made out of the scattering plane. The asymmetry Alt 1 

is defined as an electron single spin asymmetry and can, therefore, be easily obtained by 

flipping the beam helicity. While Rl and Rt are independent of photon-helicity-dependent 

phases, the interference response functions are not. As a result, the interference response 

functions can be very sensitive to phase differences generated by non-nucleonic currents and 

final state interactions. This is particularly true of Rlt 1 which can be shown to be zero in the 

PWIA. The interference response function, Rlt-, is very sensitive to the relativity included 

in the current operator, due to the various interference contributions from the charge and 

transverse current operators [T71 HE] • 
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2. Target Polarization 

It is usually assumed that electron scattering cross sections are polarized along the direc- 
tion of the three-momentum transfer q. In this case it is possible to express the cross section 
in a form where the explicit <fi dependence can be factored out of the response functions. 
In the case of a polarized target, it is often more convenient experimentally to polarize the 
target along the direction of the electron beam. The procedure for performing the appro- 
priate rotations to obtain the cross section for this situation are described in [2]. This is 
used in the calculations presented in the paper. We present asymmetries for vector and 
tensor polarization along the beam axis. The single and double asymmetries for these two 
polarizations are defined as 

v l 1Zl{T 10 , 6) + v T TZ T (Ti , 6) + v T tK-tt{Tio, 6) + v lt TI L t{T 10 , 6) 






TioS 



T _ v L TZ L (T 20 , 0) + v T TZ T (T 20 , 0) + v TT K TT (T 20 , 0) + v LT K LT (T 20} 0) 

A d — ~ 



A v 



VLT'fcLT'jTlQ, 6) + Vt'1Zt'(T 10 , 6) 

TioS 



where 



a t _ vlt'K lt ,(T 20 , 6) + v T <K T >(T 20 , 6) 
T 2 o^-> 



£ = v L n L (u, 6) + v T n T (u, 6) + v TT n T T(u, 6) + v LT n LT {u, 6) . (A2i) 



Here Ri(T w ,0) and -Rj(T 2 o,0) denote the response functions where only T w is nonzero or 
only T 2Q is nonzero. Ri(U,0) denotes the unpolarized response functions. 

3. Polarized Proton 

For unpolarized deuterons and polarized ejected protons, the proton polarization is typ- 
ically described in terms of three unit vectors 

I = ^- (A22) 

IpI 

oxl 
h = — — x- (A23) 

\q x l\ 

s = nxl. (A24) 
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Using this coordinate system allows factorization of the response functions to display explicit 
dependence on cf). However, defining asymmetries using this coordinate system produces 
asymmetries which are not properly defined for 8 P = 0° and 180°. An alternate set of basis 
vectors defined as 

V = l (A25) 

V X I 

s' = - — - (A26) 

\y x l\ 



n 



l' x s' (A27) 



eliminates this problem but results in response functions that no longer factor to give the 
explicit dependence. 
For convenience, define 

da 5 

a(S) + ha h {S) (A28) 



dE'(Kl e dVt p/ , ! s 

where u{S) is the contribution to the cross section independent of the electron helicity and 
<Jh{S) is the part of the cross section proportional to h, the single and double asymmetries 
are now defined as 

a(0) 
and 

a(0) 
where ^ = n', l',s'. 



4 = ^M (A29) 



^v = S^ (A30) 
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